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Abstract 

By applying the method of moving frames modelling one and two 
dimensional local anisotropies we construct new solutions of Einstein 
equations on pseudo-Riemannian spacetimes. The first class of solu- 
tions describes non-trivial deformations of static spherically symmetric 
black holes to locally anisotropic ones which have elliptic (in three di- 
mensions) and ellipsoidal, toroidal and elliptic and another forms of 
cylinder symmetries (in four dimensions). The second class consists 
from black holes with oscillating elliptic horizons. 
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1 Introduction 



In recent years, there has been great interest in investigation of gravitational 
models with anisotropies and applications in modern cosmology and astro- 
physics. There are possible locally anisotropic inflational and black hole like so- 
lutions of Einstein equations in the framework of so-called generalized Finsler- 
Kaluza-Klein models II and in low-energy locally anisotropic limits of (super) 



string theories [10 



In this paper we shall restrict ourselves to a more limited problem of def- 
inition of black hole solutions with local anisotropy in the framework of the 
Einstein theory (in three and four dimensions). Our purpose is to construct 
solutions of gravitational field equations by imposing symmetries differing 
in appearance from the static spherical one (which uniquely results in the 
Schwarzschild solution) and search for solutions with configurations of event 
horizons like rotation ellipsoids, torus and ellipsoidal and cylinders. We shall 
proof that there are possible elliptic oscillations in time of horizons. 

In order to simplify the procedure of solution and investigate more deeply 
the physical implications of general relativistic models with local anisotropy 
we shall transfer our analysis with respect to anholonomic frames which are 
equivalently characterized by nonlinear connection (N-connection) structures 
05 05 B is jlOj - This geometric approach is very useful for construction of met- 
rics with prescribed symmetries of horizons and definition of conditions when 
such type black hole like solutions could be selected from an integral variety 
of the Einstein field equations with a corresponding energy-momentum ten- 
sor. We argue that, in general, the symmetries of solutions are not completely 
determined by the field equations and coordinate conditions but there are also 
required some physical motivations for choosing of corresponding classes of 
systems of reference (prescribed type of local anisotropy and symmetries of 
horizons) with respect to which the 'picture' of interactions, symmetries and 
conservation laws is drawn in the simplest form. 

The paper is organized as follows: In section 2 we introduce metrics and 
anholonomic frames with local anisotropies admitting equivalent N-connection 
structures. We write down the Einstein equations with respect to such locally 
anisotropic frames. In section 3 we analyze the general properties of the system 
of gravitational field equations for an ansatz for metrics with local anisotropy. 
In section 4 we generalize the three dimensional static black hole solution to 
the case with elliptic horizon and proof that there are possible elliptic oscil- 
lations in time of locally anisotropic black holes. The section 5 is devoted to 
four dimensional locally anisotropic static solutions with rotation ellipsoidal, 
toroidal and cylindrical like horizons and consider elliptic oscillations in time. 
In the last section we make some final remarks. 
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2 Anholonomic frames and N— connections 



In this section we outline the necessary results on spacetime differential geome- 
try [|j] and anholonomic frames induced by N-connection structures f| |9|, |10| . 
We examine an ansatz for locally anisotropic (pseudo) Riemannian metrics 
with respect to coordinate bases and illustrate a substantial geometric simpli- 
fication and reduction of the number of coefficients of geometric objects and 
field equations after linear transforms to anholonomic bases defined by coeffi- 
cients of a corresponding N-connection. The Einstein equations are rewritten 
in an invariant form with respect to such locally anisotropic bases. 
Consider a class of pseudo-Riemannian metrics 

g = ga/s (u e ) du- ® duP- 

inan + m dimensional spacetime V^ n+m \ (n = 2 and m = 1, 2), with compo- 
nents 

gij + NfNfrab N*h ae 
Nfh be h ab 



9a/3 



(2.1; 



where g^ = gij (u a ) and h ao = h ab (u a ) are respectively some symmetric n x n 
and m x m dimensional matrices, iVJ = N? (u^) is a n x m matrix, and the 
n + m dimensional local coordinates are provide with general Greek indices 
and denoted u 13 = (x\y a ). The Latin indices i,j,k,... in ( |2.1| ) run values 1,2 
and a, b, c, .. run values 3, 4 and we note that both type of isotropic, x l , and 
the so-called anisotropic, y a , coordinates could be space or time like ones. We 
underline indices in order to emphasize that components are given with respect 
to a coordinate (holonomic) basis 

e„ = da = d/du^ (2.2) 

and/or its dual 

= du^. (2.3) 
The class of metrics (|2.1| ) transform into a (n x n) © (m x m) block form 

g = 9ij (u £ ) dx { ® dx> + h ab (u £ ) (5y a ) 2 <g> (5y a ) 2 (2.4) 

if one chooses a frame of basis vectors 

5 a = 5/du a = (S/dx* = di-N? (u e ) d a , db) , (2.5) 

where di = d/x % and d a = d/dy a , with the dual basis being 

S a = 5u a = (dx\ 5y a = dy a + Nf (u £ ) dx l ) . (2.6) 



The set of coefficients N = {A^ a (u 6 )} from ( |2.5| ) and ( p.6| ) could be associ- 
ated to components of a nonlinear connection (in brief, N-connection) struc- 
ture defining a local decomposition of spacetime into n isotropic directions x % 



2 



and one or two anisotropic directions y a . The global definition of N-connection 
is due to W. Barthel (the rigorous mathematical definition of N-connection 
is possible on the language of exact sequences of vector, or tangent, subbundles) 
and this concept is largely applied in Finsler geometry and its generalizations 
[||, || . It was concluded |9|, [H]] that N-connection structures are induced under 
non-trivial dynamical compactifications of higher dimensions in (super) string 
and (super) gravity theories and even in general relativity if we are dealing 
with anholonomic frames. 

A N-connection is characterized by its curvature, N-curvature, 

Q ij = d i N i ~ d i N ? + N i 9 b N ! ~ N b d b N?. (2.7) 
As a particular case we obtain a linear connection field T^ b (x l ) if N^(x l , y a ) = 

For nonvanishing values of Qfj the basis ( |2.5|) is anholonomic and satisfies 
the conditions 

S a S/3 - S/35 a = w 1 a p5 11 

where the anholonomy coefficients w^ aj3 are defined by the components of N- 
connection, 

<j = Q,w k aj = 0,w k ia = 0,w k ab = Q,w e ttb = 0, 
w a -- = -Q a - w b ■ = -d N b w b - = d N b 

We emphasize that the elongated by N-connection operators ( |2.5| ) and 
(|2.6| ) must be used, respectively, instead of local operators of partial derivation 



(|2.2|) and differentials (|2.3|) if some differential calculations are performed with 
respect to any anholonomic bases locally adapted to a fixed N-connection 
structure (in brief, we shall call such local frames as la-bases or la-frames, 
where, in brief, la- is from locally anisotropic). 

The torsion, T(5 7 ,5^) = T a ^5 ai and curvature, R (5 T , <5 7 ) 5p = R^ 1T 5 ai 
tensors of a linear connection T a p are introduced in a usual manner and, 
respectively, have the components 

T% = r% - T« ip + vffr (2.8) 

and 

-^W = ^T^ a ^y — fiyT a p S + Y v j3l Y a ipT — r^r ^ + r a p (p w' p rr . (2.9) 

The Ricci tensor is defined 

R Pl = R^ ia (2.10) 

and the scalar curvature is 

R = g^R M . (2.11) 
The Einstein equations with respect to a la-basis (|2.6|) are written 

R 

Rp-y - ~n9Pi = kTf3~ n (2.12) 
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where the energy-momentum d-tensor includes the cosmological constant 
terms and possible contributions of torsion ( |2.8|) and matter and k is the cou- 
pling constant. For a symmetric linear connection the torsion field can be 
considered as induced by the anholonomy coefficients. For dynamical torsions 
there are necessary additional field equations, see, for instance, the case of 



locally anisotropic gauge like theories [1 I 



The geometrical objects with respect to a la-bases are distinguished by 
the corresponding N-connection structure and called (in brief) d-tensors, d- 
metrics ( |2.4|) , linear d-connections and so on || |^, |10|| . 

A linear d-connection Dona spacetime V, 

D 5i 5p = T a Pl (x\y a )5 a , 

is parametrized by non-trivial horizontal (isotropic) - vertical (anisotropic), 
in brief, h-v-components, 

= (L % jki L a bk ,C l jc ,C a bc ^ . (2.13) 

Some d-connection and d-metric structures are compatible if there are satisfied 
the conditions 

D a gp~ f = 0. 

For instance, the canonical compatible d-connection 

cpa ( c r* era c fii c r~ia \ 

1 fa — \ ^ jkl ^ bki U jci U be) 

is defined by the coefficients of d-metric ( |2.4j ), gij (x\ y a ) and h ab (x\ y a ) , and 
of N-connection, Nf = iVf (x % , y b ^j , 



C T i 

L jk - 


^g m ($kg n j + Sjg nk 


era 
^ bk — 


d b N a k + l -h ac (5 k h bc 


U jc ~ 


^g ik d c g jk , 


° be — 


h ad (d c h db + d b h dc 



;2.i4) 



h dc d b Nf - h db d c Nf h 



M'lbc) ■ 

The coefficients of the canonical d-connection generalize with respect to la- 
bases the well known Cristoffel symbols. 

For a d-connection (|2.13|) we can compute the non-trivial components of 
d-torsion (|2T8| ) 

^jk kji 1 ja ^ y .jai 1 aj ja> 

be — D .bc — °bc °c6; \ z - 10 ) 

T a bt = d b N?-L%, T% = -T a H . 



rpl 

1 .jk 


rpl 

- 2 jk - 


rpi 

■ja 


= 0, T 


rpa 

1 .ij 


— -fl a - 
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In a similar manner, putting non-vanishing coefficients ( [2.13[ ) into the for- 
mula for curvature (Hip, we can compute the coefficients of d-curvature 



R (5 T , 5 y ) 5/3 — Rp'^a, 
split into h-, v-invariant components, 

R'h.jk = 5kL % hj - 8jL hk + L hj U mk - L™ k L l mj - C z ha p, a jk , 

n b.jk — °k^.bj ~ °i lj .bk^ ^.bj^.ck ~ ^.bk^.cj ~ U .£>e".jfc> 

P'lka = ®kL % jk + C l jb T b ka — (d k C l ja + L\ k C l ja — L l j k Cj a — L c ak C l jc ), 

pc a t c , /nrc rpd /o ric j_tc /-id jd /~ic jd sic \ 

r b.ka — .bk ~t~ ^ .bd 1 .ka \ u k^ .ba ~>~ ^ .dk^ .ba ^ .bk^ .da ^.ak^.bd) 

j.bc — .jb u b^.j c ~<~ ^.jb^.hc ° .jc^hbi 

Q.a o sia o , ^~ve 

D fe.cd — °d^.bc ~ °c^.bd "+" °.6c U .ed — ^.ftd^.ec- 

The components of the Ricci tensor (|2.10| ) with respect to locally adapted 
frames Q2.5Q and (|2.6|) (in this case, d-tensor) follows: 



Rij — R'ij k , Ria — 2 Pia — P~Lkai (2.16) 

p 1 p p. ft p C-c 

-T-ai -<ai r a.ibi ^ab ° a.bc' 

We point out that because, in general, 1 P ai ^ 2 P ia the Ricci d-tensor is 
non symmetric. This is a consequence of anholonomy of la-bases. 

Having defined a d-metric of type ( |2.4|) on spacetime V we can compute 
the scalar curvature ( |2.11| ) of a d-connection D, 



*R =G af3 R al3 = R + S, (2.17) 

where R = g^Rij and S = h ab S ab . 

Now, by introducing the values of ( |2.16 ) and ( |2.17| ) into equations ( [2 . 1 2| ) , 



the Einstein equations with respect to a la-basis seen to be 

Rij-UR + S^gij = feTy, (2.18) 



2 

S a b — - (R + S^j h ab = kT ab , 

i p . . 

-'cm ^ - 1 - ai j 

P — kT- 

1 la fv J- i< 
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where Ty, T a h, T ai and T ia are the components of the energy-momentum d- 
tensor field Tg 7 (which includes possible cosmo logical constants, contributions 
of anholonomy d-torsions (|2.15|) and matter) and k is the coupling constant. 
For simplicity, we omitted the upper left index c pointing that for the Einstein 
theory the Ricci d-tensor and curvature scalar should be computed by applying 
the coefficients of canonical d-connection (|2.14|) . 
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3 An ansatz for la— metrics 



Let us consider a four dimensional (in brief, 4D) spacetime V^ 2+2 ^ (with two 
isotropic plus two anisotropic local coordinates) provided with a metric ( |2.1|) 
(of signature (-,+,+,+), or (+,+,+,-), (+,+,-,+)) parametrized by a symmetric 
matrix of type 

" 9i + qi 2 h ?J + ni 2 /i 4 q x h 3 n x h 4 ~ 

g2 + q2 2 h3 + n 2 2 h 4 q 2 h 3 n 2 h A , , 

qih q 2 h 3 h 

n\h± n 2 hi h± 

with components being some functions 

9i = gi(x J ),qi = qi(x 3 ,z),rii = ni(x 3 ,z),h a = h a (x 3 ,z) 

of necessary smoothly class. With respect to a la-basis (|2.6|) this ansatz results 
in diagonal 2 x 2 h- and v-metrics for a d-metric ( |2.6| ) (for simplicity, we shall 
consider only diagonal 2D nondegenerated metrics because for such dimensions 
every symmetric matrix can be diagonalized) . 

An equivalent diagonal d-metric (|2.4j) is obtained for the associated N- 
connection with coefficients being functions on three coordinates (x l , z), 

N? = gi (x\z), Nl = q 2 {x\z), (3.2) 
N* = n^x^z), N$ = n 2 {x\z). 



For simplicity, we shall use brief denotations of partial derivatives, like a— 
da/dx l ,a! = da/dx 2 , a* = da/dz a'= d 2 a/dx 1 dx 2 , a** = d 2 a/dzdz. 

The non-trivial components of the Ricci d-tensor ( |2.16| ) ( for the ansatz 
(|3.1|) ) when R\ = R\ and S 3 = Sf, are computed 

R\ = 7^—[-(9l+92) + ^-(9 2 2 + 9W2)+^-(9'i 2 + 9192)1 (3.3) 
S\ = ^[-K + ±{hff + ±W, 0.4) 



ft h*y h$* h\ hth\ 

2 [ \h 3 J h 3 2h 4 2 2h 3 h 4 



1 _^ • 

+ 2h A [ 2h± 4 2h 3 ^ 

hi 



2h : . 



(3.6) 



The curvature scalar R ( |2.17] ) is defined by two non-trivial components 
R = 2R\ and S = 2Sf. 
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The system of Einstein equations fl2.18|) transforms into 

R\ = - K T 3 3 = -kT\, (3.7) 

S 3 3 = -kT\ = -kTI (3.8) 

= ^ 3i , (3.9) 

Pa = /tT 4i , (3.10) 



where the values of R\,S 3 ,P ai , are taken respectively from ( |3.3j ), (p^), ( |3.5| ), 

We note that we can define the N-coefficients ( |3.2| ), qi(x k , z) and rii(x k , z), 
by solving the equations ( [3.9D and (|3.10 ) if the functions hi(x k ,z) are known 



as solutions of the equations (pTB 



Let us analyze the basic properties of equations ( p.8[) -( ^T0| ) (the h-equati- 
ons will be considered for 3D and 4D in the next sections). The v-component 
of the Einstein equations ( p.7|) 



gMty-MZ)®)-****- (311) 

(here we write down the partial derivatives on z in explicit form) follows from 
( |3.4p and ( |3.8| ) and relates some first and second order partial on z derivatives 
of diagonal components h a (x l ,z) of a v-metric with a source kT\{x % ,z) = 
kT\ = /tTg in the h-subspace. We can consider as unknown the function 
h 3 {x l ,z) (or, inversely, hi(x l ,z)) for some compatible values of h4(x\z) (or 
h 3 {x l ,z)) and source Ti(x\z). 

By introducing a new variable (3 = h*Jh± the equation ( |3.11| ) transforms 
into 

1 fih* 

P * + 2^ ~ 2h ~ 2KTlh = ° (3 ' 12) 

which relates two functions j3 (x l , z) and h% (x l , z) . There are two possibilities: 
1) to define (3 (i. e. when kT\ and h 3 are prescribed and, inversely 2) to 
find h 3 for given kTi and h 4 (i. e. (3)\ in both cases one considers only "*" 
derivatives on ^-variable (coordinates x l are treated as parameters). 

1. In the first case the explicit solutions of ( |3.12| ) have to be constructed by 
using the integral varieties of the general Riccati equation M which by a 
corresponding redefinition of variables, z — > z (<j) and (3{z) — > rj (<j) (for 
simplicity, we omit here the dependencies on x % ) could be written in the 
canonical form 

where ^ vanishes for vacuum gravitational fields. In vacuum cases the 
Riccati equation reduces to a Bernoulli equation which (we can use the 
former variables) for s(z) = (3~ x transforms into a linear differential (on 
z) equation, 

f + §- 3 S - \ = 0. (3.13) 
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2. In the second (inverse) case when h 3 is to be found for some prescribed 
kTi and (3 the equation ( |3.12|) is to be treated as a Bernoulli type equa- 
tion, 

h; = -^(hf + (^ + p)h 3 (3.i4) 

which can be solved by standard methods. In the vacuum case the 
squared on h 3 term vanishes and we obtain a linear differential (on z) 
equation. 

A particular interest presents those solutions of the equation (|3.12| ) which 
via 2D conformal transforms with a factor u = u(x\z) are equivalent to a 
diagonal h-metric on x-variables, i.e. one holds the parametrization 

h-i = u(x\ z) a 3 (x l J and h 4 = u(x\ z) a 4 (V) , (3.15) 

where 03 (x l ) and a 4 (x l ) are some arbitrary functions (for instance, we can 
impose the condition that they describe some 2D soliton like or black hole so- 
lutions). In this case (3 = uj* ju) and for 7 = u~ l the equation (|3.12| ) transforms 
into 

7 7 ** = -2/€T 1 a 3 (xf) (3.16) 
with the integral variety determined by 



-4^X103(^)^17! +Ci(x i 



+ C 2 (x"), (3.17) 



where it is considered that the source Ti does not depend on z. 

Finally, we conclude that the v-metrics are defined by the integral varieties 
of corresponding Riccati and/or Bernoulli equations with respect to ^-variables 
with the h-coordinates x % treated as parameters. 



4 3D black la-holes 

Let us analyze some basic properties of 3D spacetimes V^ 2+1 ^ (we empha- 
size that in approach (2 + 1) points to a splitting into two isotropic and one 
anisotropic directions and not to usual 2D space plus one time like coordi- 
nates; in general anisotropies could be associate to both space and/or time 
like coordinates) provided with d-metrics of type 

5s 2 = g x (x k ) (dx 1 ) 2 + g 2 (x k ) (dx 2 f + h 3 {x\ z) (5zf , (4.1) 

where x k are 2D coordinates, y 3 = z is the anisotropic coordinate and 

5z = dz + Nf(x k ,z)dx\ 
The N-connection coefficients are 

N* = gi (x\z), Nl = q 2 (x i ,z). (4.2) 
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The non-trivial components of the Ricci d-tensor ( |2.16| ), for the ansatz 
( |3.1|) with /i 4 = 1 and rij = 0, R{ = R\ and P 3 j, are 



R\ = ^J-(9:+h) + ^^ + 9W 2 )+^(9[ 2 + 9m)], (4.3) 



3i 




2gi 



h** 



(4.4) 



(for 3D the component S3 = 0, 



see 



I))- 



The curvature scalar R ( |2.17|) is R = R = 2R\. 
The system of Einstein equations ( p. 18 ) transforms into 



Rl 



Si ■ 



-f2 
1 2 



(4.5) 
(4.6) 

0; the 



which is compatible for energy-momentum d-tensors with T\ 
values of R\ and P& are taken respectively from ( |4.3| ) and ( |4.4|) . 

By using the equation ( |4.6[ ) we can define the N-coefficients (|4.2|), q; t 
if the function h 3 (x k , z) and the components T 3i of the energy-momentum d- 
tensor are given. We note that the equations (|4.4|) are solved for arbitrary 



(x k , z) 



functions h 3 = h 3 (x ) and % = qi(x , z) if Y 3i = and in this case the compo- 
nent of d-metric h 3 (x k ) is not contained in the system of 3D field equations. 



4.1 Static elliptic horizons 



Let us consider a class of 3D d-metrics which local anisotropy which are similar 
to Banados-Teitelboim-Zanelli (BTZ) black holes ||T|. 
The d-metric is parametrized 



gi{x\x 2 )(d X 1 f+{dx 2 ) 2 -h 3 (x\x\t) (St) 2 , (4.7) 



5s< 



where x 1 — r / r h for r h = const, x 2 = @/ r a if r a = y |kT|| ^ and x 2 = ^ if 
Y 3 = 0, y 3 = z = t, where t is the time like coordinate. The Einstein equations 
ft4.5|) and ( [4.6|) transforms respectively into 



and 



d 2 gi 
d( X 2 ) 2 



1 <9 2 /i 3 



J_(dgi 
2gi \dx 2 



2kT1 9i = 



1 Oh* 



h 3 dz 2 \h 3 dz 
By introducing new variables 

V = g'Jgi and s = h* 3 /h 3 



-KT 3i . 



(4.9) 



(4.10) 
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where the 'prime' in this subsection denotes the partial derivative d/x 2 , the 
equations ( |4.8| ) and (|4.9|) transform into 

2 

p' + ^- + 2e = (4.11) 

and 

s* Qi = KT 3i , (4.12) 

where the vacuum case should be parametrized for e = with x % = x l and 
e = 1(— 1) for the signature 1(— 1) of the anisotropic coordinate. 

A class of solutions of 3D Einstein equations for arbitrary ^ = qi(x k , t) and 
T^i = is obtained if s = s(x 1 )- After integration of the second equation from 
( TO) , we find 

hz{x\t) = h m ( X k )exp [s {0) ( X k ) t] (4.13) 

as a general solution of the system ( J4.12 ) with vanishing right part. Static 
solutions are stipulated by qi = qi(x k ) and S(o)(x k ) — 0- 

The integral curve of ( |4.11|) , intersecting a point (xfo)^P{o)) > considered as 
a differential equation on x 2 is defined by the functions B 



p = ; 7 P ( o) 71 2 s ' e = ' ^ 

p (0) -2tanh(x 2 -xj )) n 1K . 

p = : . ^ , ; / o tt. e> °; ( 4 - 15 ) 



l + ^tanh^-^ 
P(o) " 2 tan (x 2 - X( 0) ) n 

p = — h; Hv> e <°- ( 4 - 16 ) 



Because the function p depends also parametrically on variable x 1 we must 
consider functions xfo) = X 2 ) (x 1 ) and p (0) = P( ) (x 1 ) • 

For simplicity, here we elucidate the case e < 0. The general formula for 
the nontrivial component of h-metric is to be obtained after integration on x 1 
of ( gig ) (see formula flOg)) 

01 (x\x 2 ) = 01(0) (x 1 ) jsinfx 2 - X(o) (x 1 )] + arctan p ^ x ^ j > 
for P(o) (X 1 ) 7^ 0, and 

0i (x\ X 2 ) = 01(0) (x 1 ) cos 2 [x 2 - X( ) (x 1 )] ( 4 - 17 ) 

for P(o) (x 1 ) = 0, where 5-1(0) (x 1 ) , X 2 ) (x 1 ) and p {0) (x 1 ) are some functions 
of necessary smoothness class on variable x 1 = ^/y/TiS, when e is the en- 
ergy density. If we consider = and a nontrivial diagonal components 
of energy-momentum d-tensor, = diag[0, 0, — e], the N-connection coeffi- 
cients qi{x l -,t) could be arbitrary functions. 
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For simplicity, in our further considerations we shall apply the solution 



(£T7D. 

The d-metric fl4.7|) with the coefficients ( 4.17 ) and ( 4.13|) gives a general 
description of a class of solutions with generic local anisotropy of the Einstein 



equations (|2.18|) . 

Let us construct static black la-hole solutions for S( ) (xj — in (|4.13|) . 

In order to construct an explicit la-solution we have to chose some coeffi- 
cients ^3(o)(x fc ); 01(0) (x 1 ) an d Xo (x 1 ) from some physical considerations. For 
instance, the Schwarzschild solution is selected from a general 4D metric with 
some general coefficients of static, spherical symmetry by relating the radial 
component of metric with the Newton gravitational potential. In this sec- 
tion, we construct a locally anisotropic BTZ like solution by supposing that it 
is conformally equivalent to the BTZ solution if one neglects anisotropies on 
angle 9), 



r -M + 



P- 



where M = const > and — l/Z 2 is a constant (which is to be considered 
the cosmological from the locally isotropic limit. The time-time coefficient of 
d-metric is chosen 



hs (x\ X 2 ) = r- 2 A 3 (V , X 2 ) cosV - X f 0) (V)]. 



(4.18) 



If we chose in ( 4.18 



when the constant 



-Mo + j) , 



'Mai 



defines the radius of a circular horizon, the la-solution is conformally equiv- 
alent, with the factor r _2 cos 2 [x 2 — xfo) (x 1 )]; to the BTZ solution embedded 
into a anholonomic background given by arbitrary functions qi(x l , t) which are 
defined by some initial conditions of gravitational la-background polarization. 

A more general class of la-solutions could be generated if we put, for in- 
stance, 

,2 2 



with 



h{\ l ,X 2 ) =(-M(0) +7 j) , 
M(9) M ° 



(1 + ecos^) 2 ' 

where e < 1. This solution has a horizon, A3 = 0, parametrized by an ellipse 

r h 

r = 

1 + e cos 9 



with parameter rh and eccentricity e. 
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We note that our solution with elliptic horizon was constructed for a di- 
agonal energy-momentum d-tensor with nontrivial energy density but with- 
out cosmological constant. On the other hand the BTZ solution was con- 
structed for a generic 3D cosmological constant. There is not a contradic- 
tion here because the la-solutions can be considered for a d-tensor = 
diag[pi - l/l 2 ,P2 - l/l 2 , —£ - l/l 2 } with p li2 = l/l 2 and £( e ff) = £ + l/^ 2 (for 
e = const the last expression defines the effective constant r a ). The locally 
isotropic limit to the BTZ black hole could be realized after multiplication on 
r 2 and by approximations e ~ 0, cos[# — 6q {x 1 )} — 1 and qi(x k ,t) ~ 0. 



4.2 Oscillating elliptic horizons 

The simplest way to construct 3D solutions of the Einstein equations with 
oscillating in time horizon is to consider matter states with constant nonvan- 
ishing values of T 31 = const. In this case the coefficient could depend on 
t-variable. For instance, we can chose such initial values when 

h 3 ( x \ e, t) = r - 2 (-m (t) + r ~Pj cos 2 [e - e (x 1 )] (4.19) 

with 

M = M exp (—pt) smujt, 
or, for an another type of anisotropy, 

h 3 ( X \ 6, t) = r- 2 \-M + T -Pj cos 2 6 sin 2 [£ - 6 ( X \t)] (4.20) 



with 

, , , , — , . Tn 

cosf 



1 (x 1 ; t) — e 1 y — cos wit 



when the horizon is given parametrically, 

T 

r = -coscji*, 

1 + e cos v 

where the new constants (comparing with those from the previous subsection) 
are fixed by some initial and boundary conditions as to be p > 0, and cD and 
tjj\ are treated as some real numbers. 

For a prescribed value of h^(x l , 0, t) with non-zero source T31, in the equa- 
tion (|4.6| ), we obtain 

qi { X \0,t) = kT u (^\n\h z { X \e,t)\\ . (4.21) 

A solution (|4.1|) of the Einstein equations ([4.5D and (|4.6| ) with g2(x l ) — 1 
and giix 1 ^) an d ^(x 1 ^^) given respectively by formulas ( |4.17| ) and ([4.19 ) 
describe a 3D evaporating black la-hole solution with circular oscillating in 
time horizon. An another type of solution, with elliptic oscillating in time 
horizon, could be obtained if we choose ( 4.20|) . The non-trivial coefficient of 
the N-connection must be computed following the formula ( [4.21 ). 
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5 4D la— solutions 



5.1 Basic properties 

The purpose of this section is the construction of d-metrics which are confor- 
mally equivalent to some la-deformations of black hole, torus and cylinder like 
solutions in general relativity. We shall analyze 4D d-metrics of type 

5s 2 = g x (x k ) (dx 1 ) 2 + (dx 2 ) 2 + h 3 (x\ z) (5zf + h A {x\ z) (dy 4 ) 2 . (5.1) 



The Einstein equations (3/7) with the Ricci h-tensor ( p.3|) and diagonal 



energy momentum d-tensor transforms into 

h ^ -2kT|#i = 0. (5.2) 



d 2 g l 1 / dg x x 2 



d(x 2 ) 2 2g x \dx 2 J 

By introducing a dimensionless coordinate, x 2 = x 2 / \J\kT$\, and the variable 
P = 9x1 9ii where by 'prime' in this section is considered the partial derivative 
d/x 2 , the equation ( |5.2|) transforms into 

2 

p' + y + 2e = 0, (5.3) 

where the vacuum case should be parametrized for e = with x % = x% an d 
e = 1(— 1). The equations ( |5.2| ) and ( |5.3|) are, correspondingly, equivalent to 
the equations ( (4.8| ) and ( [4.11[ ) with that difference that in this section we are 
dealing with 4D coefficients and values. The solutions for the h-metric are 
parametrized like ( f4.14|) , (|4.15|) , and fl4.16|) and the coefficient gi{x l ) is given 



by a similar to ( |4.17|) formula (for simplicity, here we elucidate the case e < 0) 



which for p( ) (x 1 ) = transforms into 

9i (x\x 2 ) = 9i(o) (x 1 ) cos 2 [x 2 - X( 0) (x 1 )]' ( 5 - 4 ) 

where g x (x 1 ) , X(o) (x 1 ) an d P(o) (x 1 ) are some functions of necessary smooth- 
ness class on variable x 1 — x 1 /^/ks^ e is the energy density. The coefficients 
9i (x\x 2 ) (|5.4| ) and ^(x 1 ,^ 2 ) — 1 define a h-metric. The next step is the 
construction of h-components of d-metrics, h a = h a (x l ', z ), for different classes 
of symmetries of anisotropies. 

The system of equations ( |3.8| ) with the vertical Ricci d-tensor component 
( |3.4|) is satisfied by arbitrary functions 

h = a 3 (x l ) and h 4 = a 4 (x l ). (5.5) 

For v-metrics depending on three coordinates (x\ z) the v-components of the 
Einstein equations transform into (|3.11|) which reduces to (|3.12|) for prescribed 
values of hs(x\ z )^ an d, inversely, to ( p.!4| ) if h^(x\ z ) is prescribed. For h- 



metrics being conformally equivalent to (|5.5|) (see transforms ( [3.15 )) we are 
dealing to equations of type ( ^.16[ ) with integral varieties ( |3.17| ) 
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5.2 Rotation Hypersurfaces Horizons 

We proof that there are static black hole and cylindrical like solutions of the 
Einstein equations with horizons being some 3D rotation hypersurfaces. The 
space components of corresponding d-metrics are conformally equivalent to 
some locally anisotropic deformations of the spherical symmetric Schwarzschild 
and cylindrical Weyl solutions. We note that for some classes of solutions the 
local anisotropy is contained in non-perturbative anholonomic structures. 



5.2.1 Rotation ellipsoid configuration 

There two types of rotation ellipsoids, elongated and flattened ones. We ex- 
amine both cases of such horizon configurations. 



Elongated rotation ellipsoid coordinates: 

An elongated rotation ellipsoid hypersurface is given by the formula |7[ 



x 2 + y 2 z 1 
n + — 



P 



(5.6) 



where a > 1 and p is similar to the radial coordinate in the spherical symmetric 
case. 

The space 3D coordinate system is defined 

x = psinhusinf cos<£>, y = psinhwsinf sin 9?, z = p cosh u cos v, 

where a = coshw, (0 < u < 00, < t> < 7r, < tp < 2n). The hypersurface 
metric is 



9uu = 9w = p 2 (sinh 2 u + sin 2 
g<p<p = p 2 sinh 2 u sin 2 v. 

Let us introduce a d-metric 

5s — gi(u, v)du + dv + h 3 (u, v, <p) (6t) + h 4 (u, v , <p) (Sip) 



(5.7) 



(5i 



where 5t and 5p are N-elongated differentials. 

As a particular solution ( |5.4| ) for the h-metric we choose the coefficient 



gi (u, v) = cos 2 v. 
The h 3 (u,v,p) = h 3 (u,v,p(u,v,p)) is considered as 

"l - 



(5.9) 



h 3 (u,v,p) 



1 



sinh u + sin v 



4p 



(5.10) 



In order to define the /14 coefficient solving the Einstein equations, for simplicity 
with a diagonal energy-momentum d-tensor for vanishing pressure we must 
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solve the equation ( |3 . 1 2| ) which transforms into a linear equation (|3.13| ) if 
Ti = 0. In our case s (u,v,ip) = (3~ l (u,v,ip) , where (3 = (dh^/dip) /h&, must 
be a solution of 



ds d\nJ\h 3 
1 

dip dip 



1 

S = 2- 



After two integrations (see 0) the general solution for h,4(u,v,<p), is 



Ii4(u, v, ip) = a>4 (u, v) exp 



— J F(u, v, z) dz 



(5-11) 



where 



F(u,v,z) = l/{yj\h 3 (u,v,z)\[s 1{0 ) (u,v) + - / y/\h 3 (u,v,z)\dz]} 



z (u,v) 



s i(o)(u,v) and zq(u,v) are some functions of necessary smooth class. We 



note that if we put = (14(11, v) the equations (|378|) are satisfied for every 

h 3 = h 3 (u,v,<p). 

Every d-metric ( |5.8| ) with coefficients of type ( |5.9| ), (|5.10|) and ( |5.11|) solves 
the Einstein equations (|3.7|) - (|3.1C1|) with the diagonal momentum d-tensor 

Tp = diag [0, 0, — e = —m , 0] , 

when r q = 2ftm ; we set the light constant c = 1. If we choose 



0,4 (u, v) 



sinh u sin 2 v 
sinh 2 u + sin 2 v 



our solution is conformally equivalent (if not considering the time-time com- 
ponent) to the hypersurface metric (|5.7| ). The condition of vanishing of the 
coefficient ( [5.10| ) parametrizes the rotation ellipsoid for the horizon 



x 2 + y 2 
a 2 -I 



z 2 



where the radial coordinate is redefined via relation r = p + . After 
multiplication on the conformal factor 



(sinh 2 u + 



sin 2 v 



1 + 



4p 



approximating g\(u, v) = cos 2 v ~ 1, in the limit of locally isotropic spherical 
symmetry, 



~2 1 ~2 1 ~2 2 

x +y + z = r 
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the d-metric 05.81) reduces to 



ds 1 



1 + ^ 
4p 



(dx 2 + dy 2 + cfc 2 ) - 



4p 



rdt 2 



which is just the Schwazschild solution with the redefined radial coordinate 
when the space component becomes conformally Euclidean. 

So, the d-metric (|5.8|) , the coefficients of N-connection being solutions of 
(|3.9|) and ( |3.10|) , describe a static 4D solution of the Einstein equations when 
instead of a spherical symmetric horizon one considers a locally anisotropic 
deformation to the hypersurface of rotation elongated ellipsoid. 



Flattened rotation ellipsoid coordinates 

In a similar fashion we can construct a static 4D black hole solution with 
the horizon parametrized by a flattened rotation ellipsoid [[7[], 

x 2 + y 2 z 2 - 
1 + o l o l 

where a > and a = sinhw. 

The space 3D special coordinate system is defined 

x = p cosh ii sin t> cosy?, y = pcoshu sint> sin 9?, z = psinhwcosf , 

where 0<m<oo, 0<f<7r, < ip < 2n. 
The hypersurface metric is 

g uu = 9w = P 2 (sinh 2 u + cos 2 vj , 
9<pip = p 2 sinh 2 u cos 2 v. 

In the rest the black hole solution is described by the same formulas as in the 
previous subsection but with respect to new canonical coordinates for flattened 
rotation ellipsoid. 



5.2.2 Cylindrical, Bipolar and Toroidal Configurations 



We consider a d-metric of type (|5.1| ). As a coefficient for h-metric we choose 
9i{x 1 iX 2 ) — ( cos X 2 ) 2 which solves the Einstein equations (p3.7|). The energy 
momentum d-tensor is chosen to be diagonal, = diag[0, 0, — e, 0] with 
e ~ m = J m/unjdl, where £ni n ) is the linear 'mass' density. The coefficient 
hz (x l , z) will be chosen in a form similar to ( |5.10| ), 



4p 



1 2 



/ 



Ap 



for a cylindrical elliptic horizon. We parametrize the second v-component 
as = a4(x\x 2 ) when the equations ( p.8|) are satisfied for every h 3 = 
h(x\x 2 ,z). 
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Cylindrical coordinates: 

Let us construct a solution of the Einstein equation with the horizon having 
the symmetry of ellipsoidal cylinder given by hypersurface formula @ 



x 2 

— + 



a 2 - 1 



P*, z = z, 



where a > 1. The 3D radial coordinate r is to be computed from p 2 = p 2 + z 2 . 
The 3D space coordinate system is defined 

x = p* cosh u cos v, y = p* sinh u sin t> sin, z = z, 



where <r = coshw, (0 < u < oo, < u < 7r). 
The hypersurface metric is 



9u 



9% 



p^ sinh it + sin v) ,g zz = 1. 



(5.12) 



A solution of the Einstein equations with singularity on an ellipse is given 



by 



1 



4p 



p 2 ^sinh 2 m + sin 2 v 
1 



1 + 

4p 



6 ■ 



04 



p 2 (sinh 2 u + sin 2 



where r — p (l + . The condition of vanishing of the time-time coefficient 
h,3 parametrizes the hypersurface equation of the horizon 



3_ r 

a 2 a 2 -l 



P*{g) 



where p* (ff) = 2Km {Hn) . 

By multiplying the d-metric on the conformal factor 



p 2 (sinh 2 u + sin^ 



1 + 



4p 



where r g — J p*( g )dl (the integration is taken along the ellipse), for p* — > 1, in 
the local isotropic limit, sinw w 0, the space component transforms into ( p-12| ). 

Bipolar coordinates: 

Let us construct 4D solutions of the Einstein equation with the horizon 
having the symmetry of the bipolar hypersurface given by the formula M 



x 2 + y 2 — 



P 



tana 



+ r 



p 2 



■ 2 ' 

sin a 
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which describes a hypersurface obtained under the rotation of the circles 



y - 



P 



tan a 



+ r 



p 2 



sin 2 a 



around the axes Oz\ because |ctana| < | sina| 1 , the circles intersect the axes 
Oz. The 3D space coordinate system is defined 



x 



p sm a cos <p 

cosh r — cos a ' 
r sinh r 



y 



p sm cr sm <p 
cosh r — cos a ' 



cosh r — cos cr 
The hypersurface metric is 



-oo < r < oo,0 < cr < 71,0 < ip < 2n) . 



9tt = 9 era 



P 2 



. 2 i Qtpf 



~9 • 2 

p sm a 



\2 ' 



(5.13) 



(cosh r — cos cr) (cosh r — cos a)' 

A solution of the Einstein equations with singularity on a circle is given by 



h. 



1-^ 
4p 



/ 



Ap 



and = a 4 = sin 2 a, 



where r = p(l + . The condition of vanishing of the time-time coefficient 
h 3 parametrizes the hypersurface equation of the horizon 



x 1 i j - - c\;\\\ a ) + z 2 



4 sin 2 a ' 



where r g = J p^ g )dl (the integration is taken along the circle), p*( s ) = 2nm^ in y 
By multiplying the d-metric on the conformal factor 



(cosh r — cos a) 



1 + % 



4p 



(5.14) 



for p* — > 1, in the local isotropic limit, sint> w 0, the space component trans- 
forms into ( [5.131 ). 



Toroidal coordinates: 

Let us consider solutions of the Einstein equations with toroidal symmetry 
of horizons. The hypersurface formula of a torus is 



x 2 + y 2 — p c tanh a) + z 



The 3D space coordinate system is defined 



2 73. 
-2 . r 



sinh a 



x 



p sinh r cos <p 

cosh r — cos a ' 
psinh cr 

cosh t — cos cr 



p sm a sm if 
^ cosh r — cos cr ' 

-7r<cr<7r,0<r<oo,0<(^< 2n) . 
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The hypersurface metric is 



fjo 



P 2 



(cosh r — cos o 



2 i 9<Pf> 



~9 • 9 

p sm a 



(cosh T — COS cr 



,2 • 



(5.15) 



This, another type of solution of the Einstein equations with singularity on 
a circle, is given by 



4p 



/ 



4p 



and /i4 = 04 = sinh a, 



where r = p (l + ^ . The condition of vanishing of the time-time coefficient 
fi3 parametrizes the hypersurface equation of the horizon 



x 2 + y 2 — 



2 tanh a 



+ r 



4 sinh a 



2 > 



where r g = j p*( g )dl (the integration is taken along the circle), = 2nm^ in y 
By multiplying the d-metric on the conformal factor (|5.14f ), for — > 1, in 
the local isotropic limit, sinv ~ 0, the space component transforms into ( p,15| ). 



5.3 A Schwarzschild like la— solution 



The d-metric of type 



is taken 



5s< 



9i(x\ mx 1 ? + dd 2 + h 3 ( X \ 0, <p) {5tf + h 4 ( X \ 0, <p) m 2 , (5.16) 



where on the horizontal subspace \ l = p/ r a is the dimensionless radial co- 
ordinate (the constant r a will be defined below), x 2 = 9 and in the vertical 



subspace y 



t and y = if. The energy-momentum d-tensor is taken to 



be diagonal Y^j = diag[0, 0, —e, 0]. The coefficient gi is chosen to be a solution 
of type ( pT4| ) 

9l ( x 1 ,e)=cos 2 e. 

For 

[1 " r a /M 2 



hi = sin 9 and hs (p) 



[1 + r a /4p] 



6 ' 



where r = p ( 1 + ,r 2 = x 2 + y 2 + z 2 , r a =r g is the Schwarzschild grav- 
itational radius, the d-metric ( 5.16|) describes a la-solution of the Einstein 
equations which is conformally equivalent, with the factor p 2 (1 + r g /Ap) 2 , to 
the Schwarzschild solution (written in coordinates (p,9,ip,t)), embedded into 
a la-background given by non-trivial values of ^(p, 9, t) and nj(p, 9, t). In the 
anisotropic case we can extend the solution for anisotropic (on angle 9) grav- 
itational polarizations of point particles masses, m = m(9) , for instance in 
elliptic form, when 

T ' a {0) = (l + ecos0) 
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induces an ellipsoidal dependence on 9 of the radial coordinate, 



P 4(l + ecos#)' 

We can also consider arbitrary solutions with r a = r a (9, t) of oscillation type, 
r a ~ sin (cJit) , or modelling the mass evaporation, r a ~ exp[— st], s = const > 
0. 

So, fixing a physical solution for h 3 (p,9,t), for instance, 

[1 -r a exp[-st]/4p(l + ecosfl)] 2 
3lP ' ' J [l+r a exp[- S t]/4p(l + ecosfl)] 6 ' 

where e = const < 1, and computing the values of qi(p,9,t) and rii(p,9,t) 
from (p.9|)/ and ( |3.10| ), corresponding to given h 3 and h^, we obtain a la- 
generalization of the Schwarzschild metric. 

We note that fixing this type of anisotropy, in the locally isotropic limit we 
obtain not just the Schwarzschild metric but a conformally transformed one, 
multiplied on the factor 1/p 2 (1 + r g /Ap) 4 . 



6 Final remarks 

We have presented new classes of three and four dimensional black hole solu- 
tions with local anisotropy which are given both with respect to a coordinate 
basis or to an anholonomic frame defined by a N-connection structure. We 
proved that for a corresponding ansatz such type of solutions can be imbedded 
into the usual (three or four dimensional) Einstein gravity. It was demon- 
strated that in general relativity there are admitted static, but anisotropic 
(with nonspheric symmetry), and elliptic oscillating in time black hole like 
configurations with horizons of events being elliptic (in three dimensions) and 
rotation ellipsoidal, elliptic cylinder, toroidal and another type of closed hy- 
persurfaces or cylinders. 

From the results obtained, it appears that the components of metrics with 
generic local anisotropy are somehow undetermined from field equations if the 
type of symmetry and a correspondence with locally isotropic limits are not 
imposed. This is the consequence of the fact that in general relativity only a 
part of components of the metric field (six from ten in four dimensions and 
three from six in three dimensions) can be treated as dynamical variables. This 
is caused by the Bianchi identities which hold on (pseudo) Riemannian spaces. 
The rest of components of metric should be defined from some symmetry 
prescriptions on the type of locally anisotropic solutions and corresponding 
anholonomic frames and, if existing, compatibility with the locally isotropic 
limits when some physically motivated coordinate and/or boundary conditions 
are enough to state and solve the Cauchy problem. 

Some of the problems discussed so far might be solved by considering the- 
ories containing non-trivial torsion fields like metric-affine and gauge gravity 
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and for so-called generalized Finsler-Kaluza-Klein models. More general solu- 
tions connected with locally anisotropic low energy limits in string/M-theory 
and supergravity could be also generated by applying the method of compu- 
tation with respect to anholonomic (super) frames adapted to a N-connection 
structure. This topic is currently under study. 
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